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ABSTRACT 

We report the detection of a 5.8 x 10 4 s periodicity in the 0.5-10 keV X-ray light curve of 
the Seyfert galaxy IRAS18325-5926, obtained from a 5-day ASCA observation. Nearly 
9 cycles of the periodic variation are seen; it shows no strong energy dependence and 
has an amplitude of about 15 per cent. Unlike most other well-studied Seyfert galaxies, 
there is no evidence for strong power-law red noise in the X-ray power spectrum of 
IRAS18325-5926. Scaling from the QPOs found in Galactic black hole candidates 
suggests that the mass of the black hole in IRAS18325-5926 is ~ 6 x 10 6 -4 x 10 7 M©. 
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1 INTRODUCTION 

ASCA X-ray spectra of many Seyfert galaxies show that 
, the iron emission line is broad and skewed to low ener- 
gies (Tanaka et al 1995; Nandra et al 1997), indicating that 
the X-rays originate from above the surface of an accretion 
disk close to a massive black hole (Fabian et al 1995). The 
line profiles demonstrate that much of the X-ray emission 
orginates from about 10-20 gravitational radii, but do not 
constrain the mass of the black hole. For that the emission 
radius is required in physical units. 

X-ray variability is an important characteristic of such 
active galactic nuclei (AGN) and may hold the key to obtain- 
ing that size. Short time scale X-ray variations are probably 
due to individual flares above the accretion disk and so do 
not determine the orbital radius unless the flare mechanism 
is understood much more than is currently the case. A pe- 
riodic signal would suffice if the origin of the period, orbital 
or otherwise, was known. 

Previous studies of X-ray variability of AGNs have 
found that power spectra are featureless and power-law in 
shape (e.g., Lawrence et al 1987; McHardy 1989). No char- 
acteristic periodicity has been found except for possible 
quasi-periodic oscillations (QPOs) in NGC5548 (~ 500s, 
Papadakis & Lawrence 1993; Tagliaferri et al 1996 argue 
against the presence of QPO in those data), NGC4051 
(~ 1 hr, Papadakis & Lawrence 1995) and RX J0437.4-4711 
(0.906 ± 0.018 dy, Halpern & Marshall 1996). Very clear 
QPO have been seen in several Galactic Black Hole Candi- 



dates; a 67 Hz QPO in GRS 1915+105 (Morgan, Remillard 
& Greiner 1997) and QPO at lower frequencies, 1-10 Hz, in 
GX339-4 and GS 1124-68 (Dotani 1992; Belloni et al 1997). 

Here we report the discovery of a periodic signal from 
IRAS18325-5926, which is a Seyfert galaxy with a broad iron 
line (Iwasawa et al 1995, 1996a). The X-ray source is one of 
the complete X-ray sample selected by the HEAO-1 A2 sur- 
vey (Piccinotti et al 1982) and it was identified with the 
IRAS galaxy (= Fairall 49) by Ward et al (1988). Although 
its narrow-emission-line-dominated optical spectrum is sim- 
ilar to that of Seyfert 2 galaxies (DeGrijp et al 1985), the 
presence of a weak broad wing to the Balmer lines (Carter et 
al 1984; Iwasawa et al 1995) suggests a dust-obscured broad- 
line region in this galaxy. The X-ray spectrum is moderately 
absorbed by a column density of ~ 1 x 10 22 cm~ 2 , which is 
consistent with the picture of an obscured Seyfert 1 nucleus. 
Apart from the slightly steep continuum slope (F ~ 2.1), the 
absorption- free hard X-ray emission (> 2 keV) has similar 
properties to those of Seyfert 1 galaxies. The X-ray source 
has been found to be highly variable, as observed in other 
Seyfert 1 galaxies such as MCG-6-30-15. 

A five-day long ASCA observation with a net exposure 
time of ~ 200 ks has been carried out with the aim of inves- 
tigating X-ray variability in the source, including the iron 
line and underlying continuum. The general spectral prop- 
erties of the X-ray source are similar to those seen before, 
and will be discussed in another paper. We report in this 
Letter on a timing analysis of the ASCA data. 
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Figure 1. The ASCA SO light curve in the 0.5-10 keV band. Each data bin is of exposure time 1000 s for clarity while the timing 
analysis used 360 s bins. The epoch of the start is 1997 March 27, 13:21:31 (UT). The nine stars plotted at the top are spaced at the 
detected frequency. 



2 THE ASCA DATA 

IRAS 18325-5926 was observed with ASCA for 5 days be- 
tween 1997 March 27 and 1997 March 31. The X-ray tele- 
scopes pointed at the source at the 1CCD mode nominal 
position. Averaged count rates detected with the four de- 
tectors over the total observing run are shown in Table 1. 
The Solid state Imaging Spectrometers (SIS; SO and SI) are 
operated with the Faint mode using one standard CCD chip 
in each detector. The standard settings of PH mode were 
used for the Gas Imaging Spectrometers (CIS; G2 and G3). 
The data reduction was performed using FTOOLS (version 
4.0) and standard calibration. After data selection, the net 
exposure times of useful data are 204 ks for the SIS and 180 
ks for the CIS. 

The source events were collected from a circular region 
centred on the source with a radius of 4 arcmin for the SIS 
and 6 arcmin for the CIS. Background in each detector was 
taken from a source-free region on the same detector. The 
backgrounds were basically stable with an intensity about 
7-8 per cent of the mean counts collected from the source 
regions. 

The mean X-ray flux of the source is 1.7 x 
10~ n crg cm~ 2 s _1 in the 2-10 keV band during the ob- 
servation, a factor of ~ 2 brighter than the previous ASCA 
observation in 1993 but a factor of ~ 1.5 fainter than the 
Ginga observation in 1989. The absorption-corrected 2-10 
keV luminosity is L2-iokeV ~ 3 x 10 43 erg s _1 (Ho = 50 
km s _1 Mpc -1 ). 

The minimum time resolution of the data is 4 s. X-ray 
light curves from the four detectors have been binned into 
360 s for each data point for the timing analysis presented 
below. The count rate decreased gradually over the whole 



Table 1. Mean count rates observed with each of the four de- 
tectors during the observation. No corrections for vignetting has 
been made. 

Detector Count rate 





cts" 1 


SO 


0.372 


SI 


0.307 


G2 


0.280 


G3 


0.332 



observation with many flares on short timescales (Fig. 1). 
The light curve appears to have a characteristic time scale 
of ~ 6 x 10 4 s as well as random shorter timescales. We 
have searched for a possible periodicity of the X-ray mod- 
ulation and examined its significance using a periodogram 
technique. 



3 TIMING ANALYSIS 

The Lomb algorithm (Lomb 1976) is appropriate for search- 
ing for a periodic modulation from unevenly sampled data 
such as inherent in observations from low orbit satellites like 
ASCA (Tanaka, Inoue & Holt 1994). In Fig. 2, the power- 
spectrum computed using the Lomb periodogram (Press et 
al 1992) is shown. 

A large peak with a period just longer than the dataset 
and due to the trend seen across the light curve, lies at 
the lowest frequency shown in Fig. 2. There are two other 
significant peaks. The higher frequency peak at 1.72 x 10~ 4 
Hz or 96.9 min is likely to originate in the orbital period of 
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Table 2. Results of the Lomb periodgram analysis. Probabilities 
that the signals could occur randomly from the four detectors 
are shown for the peak at 1.725 X 10 — 5 Hz. G2 (which had the 
lowest count rate from the source) has a more significant peak at 
a frequency 1.27 X 10~ 3 Hz, not seen in the other detectors. 

Detector Probability 
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Figure 2. The power spectrum derived from the ASCA SO light 
curve of IRAS18325-5926 in the 0.5-10 kcV band. A peak at the 
right corresponds to the ASCA orbital period. The highest peak 
at 1.725 x 10~ 5 Hz is intrinsic to the source. Large power at the 
low frequency below 3 X 10 — 6 Hz is due to the decreasing trend 
across the whole observation. Note that a frequency correspond- 
ing to the whole observation length is ~ 2 X 10~ 6 Hz. 



ASCA (Tanaka et al 1994). Another, more significant, one is 
found at 1.725 x 10 -5 Hz, corresponding to a period of 58.0 
ks. This peak is found strongly in all detectors (Table 2) 
except for G2 in which the observed count rate is the lowest 
among the four detectors thus the data are most noisy. The 
probabilities of the peak being random are shown in Table 
2. If there is significant clumping in the data sampling, then 
the derivation of the probability could be affected (Home & 
Baliunas 1986). We have verified that this is not the case 
using simulated light curves, one in which random count- 
rates (with the same mean and variance as the data) arc 
assigned to the observed times and the other in which the 
count rates were shuffled. 

The periodicity is confirmed with an epoch-folding tech- 
nique. The 90 per cent statistical error to the 58.0 ks peak 
for the observed data in the x 2_ (trial period) diagram, when 
it is fitted with a gaussian, is ~ 0.3 ks. The systematic error 
in deriving a coherent periodicity is found to be ~ 0.2 ks, 
applying the same technique to simulated sine waves with a 
period of 58.0 ks, holding the same time locations and sim- 

© 0000 RAS, MNRAS 000, 000-000 



CD 

b 




6 



i , , , , 1 , , , , 1 ,_i 

5xl0 4 10 5 

Time (s) 

Figure 3. The ASCA SO light curve folded at the detected period. 
Two cycles are shown for clarity. The original data (dotted line) 
are averaged into 10 bins per cycle (thick solid line). Error bars 
to count rates are calculated from the square root of the number 
of counts integrated over each bin while the standard deviation 
in each bin is typically ~ 0.12 cts -1 . 



ilar mean count-rate and amplitude to the real observation. 
Thus the error of the period is estimated to be ~ 0.4 ks. 

The folded light curve at the detected period is shown 
in Fig. 3. The modulation is arch or sinusoidal in shape with 
an amplitude of ~ 15 per cent. There is no strong evidence 
for X-ray colour variation associated with the periodicity. 

We note that the power spectra of the variability of 
Seyfert galaxies is typically power-law in shape over the fre- 
quency intervals shown in Fig. 2 (McHardy 1989). The pos- 
sibility of such red noise complicates the true assignment of 
significance to the peak in Fig. 2. The probabilities in Table 2 
are strictly for a sinusoidal signal in white noise. Apart from 
a large peak at the lowest frequencies due to the trend in 
the lightcurve, we do however see no evidence for red noise. 
Over the frequency range 3 x 10~ 6 Hz to above 10 -4 Hz, the 
peaks in the power spectrum are similar to those expected 
from white noise, apart from the spike at 1.73 x 10 -5 Hz. If it 
is argued that the spike is due to red noise, then we must ask 
why there are no similar large peaks at 10~ J Hz and below. 
We find that the power spectrum (and lightcurve) are not 
at all similar in shape to simulations made with red noise. 
(The light curve of MCG-6-30-15 from the previous ASCA 
long look (~ 4 dy, Yaqoob et al 1997) does show red-noise.) 

We made 1000 simulated lightcurves with the observed 
time sequence and (phase) randomized f~ 1 noise, normal- 
ized to give the same variance as the real dataset. The min- 
imum power (in the same units as Fig. 2) in the frequency 
range from 3 x 10~ 6 Hz to 1.5 x 10" 5 Hz was 42, decreas- 
ing to 25 when the lower frequency was 4 x 10 -6 Hz. If the 
power varies as / _1 ' 5 , which is more typical of AGN, then 
these values rise to 80 and 60 respectively. All of these far 
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Figure 4. The ROSAT PSPC 0.5-2 keV contour image overlaid 
on the digitized UK Schmidt survey plate of the IRAS18325-5926 
region. 



exceed the power seen in that frequency range in the real 
dataset. If the f~ x power is scaled down to match the ob- 
servations then the maximum power seen in the simulation 
beyond 1.5 x 1CT 5 is less than 10. We conclude that the ob- 
servations are inconsistent with any simple power-law red 
noise model. 

The peak at 1.725 x 1CT 5 Hz is nearly 10 times the fre- 
quency of the peak due to the ASCA orbital period. We 
consider this to be a chance coincidence. To the best of 
our knowledge, there is no corresponding frequency origi- 
nating in the spacecraft. This has been checked using the 
light curves of the calibration source of the two GIS detec- 
tors during the present observation and the light curve of 
two-day long ASCA observation of a stable X-ray source, 
the Centaurus cluster, kindly provided by S. Allen. An ex- 
amination of the X-ray image rules out a possibility that 
wobbles of the spacecraft causes the observed X-ray modu- 
lation. 

Note that the X-ray source is the brightest both in 
the ASCA and ROSAT PSPC images of the region and is 
located at the optical position of the galaxy (18 h 36 m 58 s , 
-59°24'09", J2000). Simultaneous ROSAT and ASCA ob- 
servations of the source in 1993 (Iwasawa et al 1996) show 
a correlated decrease in flux by a factor of about two, con- 
firming that the ROSAT and ASCA sources are the same. 
We show in Fig. 4, the ROSAT PSPC image overlaid on the 
digitized UKSTU optical image. The offset of the X-ray and 
optical position is ~ 10 arcsec, which is within the PSPC 
positional uncertainty. 

To eliminate further any residual doubt that the vari- 
ability reported here is from IRAS18325-5926, since the pe- 
riod is in the range of that of low-mass X-ray binaries and 
some cataclysmic variables, we have estimated the probabil- 
ity of a chance alignment with a Galactic object by scaling 



from the Piccinotti et al (1982) X-ray sample. This implies 
about 100 Galactic objects in the Sky above |6| = 20° at 
the flux level of IRAS18325-5926 and thus probabilities of 
about 3 x 10~ 6 and 10~ 7 of a further object in the 1 arcmin 
and 10 arcsec radius error boxes obtained from the ASCA 
and ROSAT images, respectively. We note that both images 
are consistent with a single point source. As mentioned in 
the Introduction, the observed X-ray absorption is fully con- 
sistent with the optical classification of the Seyfert galaxy, 
but unlikely for a Galactic source. We conclude that it is 
most unlikely that we are repeating the case of NGC6814, 
where the X-ray variability detected by non-imaging detec- 
tors was due to a cataclysmic variable 40 arcmin from the 
galaxy (Madejski et al 1993). The IRAS galaxy is by far the 
most probable X-ray source. 



4 DISCUSSION 

The periodicity found in the light curve from the long ob- 
servation is about 58 ks which is probably intrinsic to the 
source whilst the 97 min peak is due to the ASCA orbital 
period. A similar X-ray flux variation is apparent in the pre- 
vious ASCA PV (Iwasawa et al 1996a) and Ginga (Iwasawa 
et al 1995; Smith & Done 1996) observations. The duration 
of each of those observations was about one day and so only 
slightly longer than the period (~16 hr). The data sets are 
too widely spaced to phase together. 

The (Newtonian) orbital period at 5Rs = 10r g is 3.4 x 
W 4 MgRi s where the mass is in 10 8 MsMq and the radius 
is 10-Rirg. So it could easily be occurring at 10r g for a 1.7 x 
10 8 M© black hole or at 20r g for a 2.1 x 10 7 M© one. In other 
words, a 16 hr period is reasonable for orbital variations at 
10-20 gravitational radii in a disk around a black hole of 
mass 2 x 10 7 — 10 s M©. The Eddington limit for such masses 
is above 10 45 erg s _1 , well above that observed. 

What we do not know from the 9 cycles that we have 
seen is whether the variation is strictly- or quasi-periodic. 
Further observations are necessary to test whether the ob- 
served periodicity is long-term or part of a QPO. The qual- 
ity factor of the current variation is at least 30 and could 
be much higher. A strict periodicity will be difficult to ex- 
plain. The most obvious clock would be an object in orbit 
about the accreting black hole. The lifetime against spiral- 
in due to gravitational radiation losses would be only weeks 
for another black hole of similar mass to the primary (say 
10 7 M©), increasing as the inverse mass for much smaller ob- 
jects. A star orbiting within an accretion disk does isolate 
a particular radius (see discussion by Syer, Clarke & Rees 
1991) and could last 100,000 yr but it is not obvious how it 
could modulate the X-ray emission. 

A transient periodicity could be due to a very long-lived 
flare, or flare patch, on the disk. The appearance of small 
variation originating at a few Schwartzschild radii may be 
boosted relativistically by orbital motion (see e.g., Boiler et 
al 1997). 

A 67 Hz QPO has recently been discovered with 
the Rossi X-ray Timing Expolorer (RXTE) in the Galac- 
tic black hole candidate GRS 1915+105 (Morgan, Remil- 
lard & Greiner 1997). If this and the periodicity seen in 
IRAS18325-5926 are related to the orbital period at the 
same radius (in gravitational units) in each system, then the 
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masses of their black holes should scale as 67/1.73 x 10 5 
which is about 4 million. Assuming that the black hole in 
GRS 1915+105 is 1OM then we find that the mass of the 
black hole in IRAS18325-5926 is ~ 4 x 10 7 M©. We note that 
both GRS 1915+105 and IRAS18325-5926 are observed at 
moderately high inclinations, of about 70 deg and 50 deg, 
respectively. 

QPO at lower frequencies, 1-10 Hz, have been seen in 
the Galactic Black Hole candidates GX339-4 and GS 1124- 
68 (Dotani 1992; Belloni et al 1997). The origin is again 
unknown, but simple scaling would then indicate a mass of 
4 x 1O 6 M0, dependent on the black hole mass of the Galac- 
tic sources. In these cases, the QPO are seen mostly when 
the source is very luminous, possibly close to the Edding- 
ton limit, although they were also seen in GS 1124-68 four 
months after its peak. 

It is unclear why any particular radius is favoured. One 
possibility is the disk oscillation model proposed by Nowak 
ct al (1997) for GRS 1915+10. General relativistic effects 
cause the epicyclic frequency for matter in the disk to have 
a maximum close to, but not at, the inner radius of the disk. 
This can lead to g-mode oscillations in the disk causing QPO 
at that frequency. However, as discussed by Nowak et al 
(1997) the power in such QPO should only be a few per cent 
of that of the disk emission. It is unlikely that they cause a 
10 per cent oscillation of the emission from the corona above 
the disk. 

Perhaps the spin of the black hole is responsible, if its 
axis is offset from that of the disk. The Bardeen-Petterson 
effect (1975) due to the dragging of inertial frames close 
to the black hole then causes the inclination of the disk to 
change at a radius of some tens of gravitational radii (see 
Rees 1984). If the inclination of the disk is moderately high 
then oscillations near that radius may be detectable through 
obscuration and/or reflection. 
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